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Introduction 

One  model  for  studying  the  dynamics  of  a  contained 
plasma  is  called  ideal  magnetohydrodynamics .   The  plasma  is 
assumed  to  be  described  by  pressure,  velocity,  and  density. 
The  plasma  is  assumed  to  be  in  local  thermodynamic  equilib- 
rium, and  to  be  governed  by  the  laws  of  conservation  of  mass, 
energy,  and  momentum.   Lorentz  forces,  J  x  B,  are  present 
but  no  electric  forces  are  included.   Forces  arising  from 
the  electric  field  are  ignored  since  density  is  assumed  to 
be  zero  in  a  plasma.   Finally,  the  plasma  is  assumed  to  be 

perfectly  conducting,  i.e.,  E  +  u  x  b  =  0,  so  that  B  is  de— 

%   'V   'b  a. 

termined  by  the  Faraday  law,  B,  =  -V  x  e  =  Vx(u  x  b) .   In 

%t     ^        %        ^     %        % 

this  model,  displacement  current  is  ignored  and  J  =  V  x  b. 

'\,    %   '\^ 

In  this  paper,  we  study  the  equilibrium  and  stability 
of  a  system  consisting  of  a  plasma  column  of  length  H,  ap- 
proximately a  cylinder  of  radius  r  ,  surrounded  by  a  vacuum 
region,  all  enclosed  within  an  outer  wall  approximately  a 
perfect  conducting  cylinder  of  radius,  pr  .   Vacuum  magnetic 
fields  B"'"^  and  B°^  are  present  both  in  the  plasma  and 
vacuum  regions,  while  surface  currents  are  present  on  the 
plasma  vacuum  interface. 

In  equilibrium,  the  stress  on  the  plasma  vacuum  inter- 
face from  the  plasma  equals  the  stress  from  the  external 
magnetic  field.   In  this  study,  the  plasma  fluid  pressure,  P; 


Is  constant,  characterized  by  a  dimenslonless  parameter,  B, 

referred  to  as  the  plasma  beta  and  defined  as  the 

ratio   of   the   plasma  fluid  pressure   to   the   total 

pressure  both   fluid   and  magnetic.    Furthermore, 

we   assume   that   the   limiting   values   of   B    and 

B°        on  the  plasma  vacuum  Interface  are  tangential  to  the 

surface,  i.e.  n  •  B  =  0.   Thus  the  condition  for  equilibrium 
% 

is  [p  +  ^  B  ]  =  0,  where  [   ]  denotes  the  magnitude  of  the 
discontinuity  of  a  function  across  a  surface.   This  equilib- 
rium is  called  a  free  boundary  equilibrium. 

We  study  the  stability  of  this  equilibrium  by  following 

the  analysis  of  Bernstein,  et  al.  (Ref.  1),  modified  by  Lust, 

et  al.  (Ref.  2).  The  dynamic  equations  of  the  plasma  system 

are  linearized  by  expressing  the  velocity,  pressure,  and 

magnetic  field  as  the  sum  of  their  values  at  equilibrium, 

u  =  0,  B  ,  p  ,  and  their  perturbed  values,  or  first  vari- 
*^o     '  %o'  '^o'  ^ 

ations,  u-,  ,  B,  ,  p,  ,  retaining  only  those  terms  linear  in  the 

perturbations.   When  the  Lagrangian  displacement,  ^,  is 

introduced,  where  u  =  5, ,  the  linearized  equations  can  be 

rewritten  as  a  formally  self-adjoint  linear  operator,  L, 

mapping  E,   into  its  second  derivative,  i.e.  E,        =   LC-   If  this 

equation  is  multiplied  by  E,.    and  Integrated  over  the  entire 

^t 

system,  the  result  is  an  expression  for  conservation  of 
energy.   The  second  variation  of  the  potential  energy,  6W, 
is  expressed  here  by  the  quadratic  form-(C,  L^) • 


The  system  is  called  stable  if  C  does  not  grow  exponen— 
tially  in  time  in  an  L^  sense.   If  5W  >  0  for  all  E, ,    then 
the  system  is  stable  in  this  sense.   The  system  is  called 
unstable  if  there  exist  C  that  grow  exponentially  in  an  L„ 
sense,  which  is  the  case  if  there  exist  ^  for  which  6W  <  0 . 

In  analogy  with  the  theory  of  linear  operators,  we 
expect  stationary  states  of  6W  to  occur  for  eigenvalues  of 
L.   Thus,  in  keeping  with  the  above  definition,  we  will  con- 
sider the  system  stable  if  all  eigenvalues  are  non^iegative , 
whereas  the  system  is  considered  unstable  if  there  exist 
negative  eigenvalues. 

In  this  study,  we  will  show  that  for  a  certain  class  of 

equilibria  characterized  by  p  and  other  parameters,  there 

exists  a  value  of  6  referred  to  as  beta  critical,  3  ,  such 

c 

that  for  all  6  <_  B  ,  the  system  is  stable,  whereas  for 


e  ,  there  exists  a  state,  K         ,    for  which  6W  =  0.   This  state 
is  referred  to  as  marginally  stable. 

In  this  study,  we  will  show  that  there  exists  a  new 
family  of  equilibria  in  a  neighborhood  of  X        •      This  is  an 
extension  to  plasma  systems  of  the  familiar  concept  of  bifur- 
cation encountered  in  certain  systems  governed  by  non— linear 
differential  equations. 


We  will  describe  our  equilibrium  as  a  perturbation  of 
system  called  a  theta  pinch.   In  a  theta  pinch,  the  magnetic 
field  is  directed  along  the  axis  of  the  column.   In  this 
study,  we  utilize  predominantly  theta  pinch  fields,  to  which 
we  add  a  small  admixture  of  fields  which  depend  on  the  usual 
cylindrical  coordinates  (|)  and  z  only  through  the  combination 
9  =  (f)  -  kz.   Such  fields  are  called  helically  symmetric,  and 
k  is  called  the  helical  wavenumber.   These  fields  will  be 
expanded  as  a  Fourier  series  of  terms  of  the  form 
r   exp  [ILG]  where  L  is  an  integer  and  m  >_  L  -  1  to  insure 
regularity.   For  our  study  we  consider  fields  consisting  of 
a  small  quantity  of  L  =  2  component  whose  ratio  to  the  theta 
pinch  field  is  characterized  by  a  small  parameter  c,  and  an 

even  smaller  amount  of  L  =  4  field,  whose  ratio  to  the  theta 

2 
pinch  field  is  characterized  by  Xc    . 

We  assume  long  helical  wavelengths,  and  introduce  a 

second  small  dimenslonless  parameter,  e,  defined  as  the 

plasma  column  radius  times  the  helical  wavenumber,  i.e., 

e  =  kr  .   We  are  then  able  to  facilitate  the  mathematical 
o 

analysis  by  expanding  all  geometrical  and  physical  quantities 
as  power  series  in  e   and  c,  and  carrying  out  a  perturbation 
analysis  of  these  expansions,  as  first  c   and  then  c  is  made 


small.   First,  e  is  made  small  enough  to  neglect  all  terms 

2 

smaller  than  order  e  .   In  the  subsequent  analysis,  e  is 

considered  fixed,  leaving  c  as  the  remaining  arbitrary,  small 


parameter.   Then  c  is  made  small  enough  to  retain  only  those 
leading  order  terms  of  significance  to  the  analysis  which 
follows . 

Again,  all  terms  past  the  leading  order  terms  are 
assumed  to  be  negligible.   We  emphasize  that  although  we 
have  not  attempted  to  justify  this  assumption,  our  results 
depend  strongly  on  its  validity. 

As  mentioned  above,  the  plasma  pressure  is  characterized 
by  the  plasma  beta.   If  the  magnetic  field  in  the  plasma 
region  were  constant,  since  the  fluid  pressure  is  constant, 
the  pressure  would  be  effectively  defined  by  specifying 
3.   However,  in  our  system,  the  magnetic  field  is  only 
approximately  constant,  and  correspondingly  we  define  an 
approximate,  constant  B,   which  merely  characterizes  rather 

than  defines  the  plasma  pressure.   For  our  study,  6  will  be 

2 
considered  small,  of  order  c  ,  appropriate  to  an  experimental 

device  called  a  Stellarator. 

The  displacement,  Cj  is  expanded  in  a  Fourier  series  of 
terms  of  the  form  c     ^r"^  exp  [KnG  -  £z)],  where  n  and  I 
are  the  azimuthal  and  axial  wavenumbers  respectively,  n  takes 
on  only  Integer  values,  m  >_  n  -  1  to  insure  regularity,  and 
the  coefficients  c     .    are  expressed  as  power  series  in  e 
and  c,  as  mentioned  above. 

For  any  value  of  n,  and  for  all  small  betas,  we  will 
show  that  there  exists  a  family  of  unstable  eigenstates. 


These  instabilities  were  also  observed  by  Freidberg  (Ref.  3) 
at  large  beta,   and  are  called  "kink  instabilities."   Since 
large  values  of  n  have  not  been  observed  in  plasma  systems, 
we  will  show  that  kink  instabilities  can  be  omitted  by 
assuming  n  <_  N  and  picking  H  correctly. 

Once  these  instabilities  are  omitted,  we  will  show  that 
there  exists  a  critical  value  of  beta,  as  described  above, 
and  correspondingly  a  marginally  stable  state  which  serves 
as  the  starting  point  for  bifurcation. 

Bifurcation  is  a  familiar  concept  encountered  in  certain 
systems  governed  by  non-linear  differential  equations  such 
as  occurs  in  the  buckling  of  beams.   When  the  physical  param- 
eters of  the  system  of  a  beam  under  compression  reach  a  state 
called  the  bifurcation  point,  the  equilibrium  representing 
the  straight  beam  becomes  unstable  to  any  further  compression 
while  another  family  of  equilibria,  representing  the  buckled 
state  of  the  beam,  becomes  the  stable  configuration.    As  the 
beam  is  compressed  further,  it  continuously  branches  to  the 
more  stable,  buckled,  configuration.   It  has  been  proven 
(Ref.  4),   in  the  theory  of  differential  equations,  that  for 
a  conservative  system,  derivable  from  a  Lagrangian,  the  bi- 
furcation point  is  always  a  marginally  stable  state  of  the 
bifurcating  equilibrium. 

More  recently,  bifurcation  analyses  of  hydromagnetic 
systems  were  studied,  first  by  Friedrichs  (Ref.   5),  and  then 


by  Yeh  (Ref.   6).   Friedrichs  studied  the  stability  of  heli- 
cal equilibria  with  small  amplitude  corrugations,  to  distur- 
bances in  the  limit  of  high  axial  and  azlmuthal  wavenumbers 
and  found  these  equilibria  unstable.   Yeh  also  studied  a 
helical  system  with  small  amplitude  corrugations,  but  found 
stable  bifurcated  equilibria  for  n  =  1  azlmuthal  wavenumbers 

In  this  study,  we  demonstrate  the  existence  of  bifur- 
cated equilibria  that  are  stable  at  states  for  which  the 
original  equilibria  are  unstable,  provided  the  outer  con- 
ducting wall  is  placed  within  1.15  r  of  the  plasma  vacuum 
interface.   We  also  show  that  the  L  =  ^  fields  do  not  affect 
the  bifurcated  stability. 


Equilibrium 

In  this  section  we  obtain  a  family  of  equilibria  con- 
sisting of  approximately  L  =  2  helically  symmetric  vacuum 
magnetic  fields  with  small  admixture  of  L  =  4  fields  for  the 
outer  region,  a  small  plasma  beta,  with  plasma  vacuum  Inter- 
face slightly  distorted  from  a  cylinder.   At  the  end  of  this 
section  we  append  some  geometrical  preliminaries  useful  in 
the  remainder  of  the  thesis.   This  section  and  the  next  one 
provide  no  new  results,  but  is   basis  for  the  bifurcation 
treatment  to  follow.   The  equilibrium  will  be  expanded  In 
terms  of  the  small  parameters,  e  and  c.  Introduced  in  the 

previous  section. 

We  first  formulate  the  helical  equilibrium  in  a  twisted 

coordinate  system  where  r,  0,  and  oi  are  unit  vectors  in  the 

directions  Vr,  V0  =  VA  -  kVz,  and  r  x  §   0  =  (}>  -  kz  being 

the  helical  angle  Introduced  in  the  previous  section.   Thus 

r,  0,  and  w  are  related  to  the  cylindrical  coordinates 

r,  (J),  and  z  by       ^        ^ 


r 

= 

r 

0 

= 

i(i- 

krz) 

03 

= 

§(kr$ 

+    z) 

2         2  2 

where  a   =  1  +  k  r  . 

Since,  from  Maxwell's  equations,  the  fields  are  diver- 
gence free,  we  can  express  the  helically  symmetric  vacuum 

field  in  terms  of  a  flux  function,  ijj(r,0),  such  that 

B-V(|;  =  0  (1) 

and  B   =    Vijjx^+Bcj  ,  (2) 

%        a,  a  w 

where    in   view    of    Maxwell's    equation,    J    =    V    x   B   =    0 

'\.        %        a, 

B      =   ^ 
oj        a 


for  C  constant,  and  ijj  satisfies 

Explicitly,  (2)  becomes 

B   =  i  ^ 

B   =  -  i^ 
^e     a  3r 

B  =  i 
w   a 

For  reasons  described  in  the  previous  section,  we  wish 

to  expand  B  in  the  small  parameter  helical  wavenumber,  k.   In 
a, 

order  to  proceed  non-dimensionally ,  we  introduce  the  small 
parameter,  e,  defined  as 


specified  precisely  later. 

Following  an  idea  of  Weitzner  (Ref.   7 ) ,  we  scale  r  and 
4;  by  setting  r  =  rr  and  ijj  =  er  i|j   and  taking  c  to  be 
independent  of  e  so  that 

Br'?.')  =t!f-^ 

B.(?,6)  =  -e  ii-  +  0(e3) 
o  ar 

B,^  =  c(l  -  e^  |-)  +  0(eS 
Equivalently ,  in  cylindrical  coordinates 

I 

B,(?,0)  =  £(C?  -  ^)  +  0(e3)  (4) 

<t)  9r 


!  (?,e)  =  C  +  e^(?  1^  -  C?^  +   (eS 


while  (3)  transforms  to 

Z2 


r  8r    9r     -2  ^q2 


We  may  express  the  above  relations  In  two  different 
ways  appropriate  for  different  purposes.   First,  If  we  work 

In  the  r,9  plane  and  define  the  usual  operators  V  and  V*  In 

_  t 

that  plane,  and  the  vector  r  In  that  plane,  then  i) 

satisfies 


V-V  i|;'  =  2C  +  OCe^)  (  5) 


and  from  ( 4 ) 


Second,  we  may  work  In  the  ?,({), z  space  with  the  conventions 
that  only  the  variables  r  and  4)  -  kz  appear  and  that  V  and 
^'    are  the  usual  derivative  operators  In  the  r  and  (}>  plane 
only,  and  r  Is  the  position  vector  In  this  plane  only.  In 
which  case  {k)    and  (5)  hold,  and  the  magnetic  field, 
expressed  by  (4),  has  the  concise  representation 

B  =  z[c  +  e^(i'-Vii;-Cr-r)]+0(e^)-ez  x  [V(4j'-  \   r-r)+0(£^)]   (7) 

In  (7) J  we  may  extend  V  to  the  derivative  In  the  full 
r,(f),z  space,  provided  r  Is  still  the  position  vector  In  the 
r,(|)  plane.   While  the  various  representations  seem  peculiar 
at  present,  they  are  particularly  appropriate  to  facilitate 
the  various  transformations  of  coordinates  which  will  be 
described  In  detail  at  the  end  of  this  section. 

A  free  boundary  problem  Is  usually  solved  by  giving  the 
external  field  sources  and  determining  the  plasma-vacuum 
interface  and  fields.   However,  here  we  employ  an  inverse 


10 


procedure  In  which  we  specify  a  flux  function  for  the  plasma 
region  satisfying  (5)-   The  flux  function  for  the  vacuum 
region  is  then  determined,  satisfying  (5)  and  the  boundary 
conditions  for  free  boundary  equilibria,  namely,  that  the 
flux  function  be  constant  at  the  plasma  vacuum  interface  as 
required  by  (1),  and  be  chosen  so  as  to  achieve  a  pressure 
balance  across  the  interface. 

This  equation  and  boundary  conditions  constitutes  a 
Cauchy  problem  for  which  the  Cauchy-Kowalewski  theorem  guar- 
antees  a  solution  for  ij;,  the  plasma  vacuum  interface,  and 
the  outer  flux  surface  in  some  neighborhood  of  the  plasma 
vacuum  interface  as  e  approaches  zero  for  small  but  finite  c , 
In  general,  a  singularity  in  the  solution  appears  within  the 
vacuum  region  for  finite  c.   However,  as  c  approaches  zero, 
the  singularity  moves  beyond  the  outer  flux  surface.   We  are 
then  able  to  justify  our  expansion  techniques  by  employing  a 
known  existence  theorem  for  free  boundary  equilibria 
(Ref.   8).   Once  we  have  solved  the  inverse  problem,  this 
theorem  can  be  used  to  solve  the  direct  problem  in  which  the 
outer  flux  surface  is  taken  as  the  initial  surface  for  the 
exact  solution  of  the  free  boundary  equilibrium  and  the 
plasma  vacuum  interface  determined.   It  is  expected  that 
this  interface  will  be  close  to  the  one  computed  for  the 
inverse  problem. 

We  proceed  to  construct  the  equilibrium  using  this 
inverse  procedure. 
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We  take  the  flux  function  for  the  plasma  region  to  be 
_2       _2  _^ 

(Jj^"^  =  y(  I ec  |-  cos  29  +  \ec^    r   cos  ^19)  (8) 

which  Is  a  solution  of  (5)  with  the  constant  y  replacing  ?. 

For  the  subsequent  analysis,  41  will  be  written  without 
primes . 

From  (7),  we  determine  the  field  In  the  plasma: 

B^"^  =  y[z  +  ecr(sln  29r  +  cos  29$) 

?  -3  -  (9) 

-  4Xec   r  (sin  49r  +  cos  49(())] 

The  radius  of  the  plasma-vacuum  Interface  is  expanded 

In  a  Fourier  series: 

r  =  1  +  a2  cos  29  +  a^^  cos  49  +  ag  cos  69  +  ...      (10) 

From  the  condition  that  ij;''"'^  be  constant  at  the  plasma 

surface,  we  solve  for  a2,  a^,  and  ag  order  by  order  In  c. 

The  result  is 

-2  =  1^  ^It-  !  ^)^^ 

a^^  =  (^  -  X)c^  (11) 

For  the  vacuum  region,  since  \p  must  also  satisfy  (5), 

we  can  with  all  generality,  set 

_2  2  _-2 

ijj^^   =  C(?—  +ec^  log  r  -  Cp  £cr   cos  29  +  c_2  ecr    cos  29 

p   ^  p  _-4 

+  c^  ec   r   cos  49  +  c_^    ec   r   cos  49  +  ...)   (12) 

We  take  c^  =  0  so  that  the  net  current  in  the  ^ree 
boundary  is  zero.  The  inclusion  of  c^  is  possible,  but 
experience  suggests  that  it  leads  to  less  stable  equilibria. 


12 


Since  the  plasma  pressure  supported  by  the  field,  at 
equilibrium,  is  just  [  p—  ],  the  plasma  beta  is 

a    -  P      _   [B^] 


p  +  B^  ^"    b2  °^t 


and  since,  from  (6) 


[b2]  =  ^2  _  ^2  ^  Q(^2^^ 


6  =  (1  -  i^)'^)    (1  +  Oie^)    . 
Without  serious  error,  we  define 

3=1-  (^)^   .  (13) 

The  remaining  boundary  conditions  for  free  boundary 
equilibria  are, 

,OUt/-v         ,    , 

\li        (r)  =  constant 


and. 


2 

[B  ]  =  constant 


at  the  plasma  surface. 

These  conditions,  along  with  (6),  (10),  (11),  (12),  and 

(13)  and  the  assumption  that  plasma  beta,  3,  is  small  of 

2 
order  c  ,  enable  us  to  solve  order  by  order  in  c  for  the 

coefficients  Cp,  c  „,  c^y    and  c_2,  of  (12).   After  a  long 

computation 

=  £   c^ 
^2    2    ~   T~ 


c6 


Cj,  =  \c      +   0(c) 


0(cS 


(1^) 
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Replacing  these  coefficients  in  (7),  we  have  for  the 
outer  vacuum  fields, 
B°^^  =  c  (z+cer(sln  2er+cos  26$)-  |^  er(sin  2er+cos  29^) 

+  e^l  f   (sin  2er  -  cos  20$)-4Xc^  r  (sin  4er+cos  4ei) 

+  (higher  order  terms  In  e   and  c)    .  (15) 

Thus,  the  inverse  procedure  has  yielded  the  desired 
vacuum  magnetic  fields. 

In  addition,  we  see  from  (10)  that  the  equilibrium  con- 
figuration of  the  plasma  surface  is  that  of  a  cylinder  with 
small  helical  corrugations  of  width,  order  c,  with  respect 
to  the  characteristic  radius,  r  ,  of  the  cylinder. 

We  conclude  this  section  by  giving  some  useful  geomet- 
rical preliminaries,  associated  with  the  several  coordinate 
transformations  Introduced  earlier  in  this  section. 

We  define  rectangular  coordinates  corresponding  to  the 
r,  9  coordinates  defined  in  this  section,  as  follows: 

X   =  r  cos  9  =  —  cos  kz  +  ^  sin  kz 
o  o 

y   =  r  sin  9  =  —  sin  kz  +  ^  cos  kz 
^o  ^o 

and  extend  these  coordinates  from  the  x  ,  y  -plane  to  a 

t   I   I 

transformation  of  x,y,z-space  to  x  ,y  ,z  -space  as  follows: 

r  X   =  X  cos  kz  +  y  sin  kz 

r  y   =  -X  sin  kz  +  y  cos  kz  (l6 ) 
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These  relations  emphasize  that,  by  studying  long 
helical  wavelengths,  we  are  assuming  slow  z  dependence 
relative  to  x  or  y  dependence.   If  we  define  the  corre— 
spending  system  of  orthogonal  vectors  x  ,  y  ,  z   as 

X   =  X  cos  kz  +  y  sin  kz 

y   =  -X  sin  kz  +  y  cos  kz  ,       . 

z   =  z 

I         I  t 

the  X   and  y   vectors  are  not  Independent  of  z  ,  In  fact 

9x'    -'  ,       8y'      ^' 

— r  =  y       and      -S-  =  -  x 
8z  8z 

so  that  when  we  Introduce  the  standard  polar  coordinates, 

II  ft 

r  ,  9  ,  In  the  x  ,  y   plane,  previously  referred  to  as 

r,  6,  and  the  unit  polar  vectors  r  ,  6  ,  we  see  that 

^=  ^=  0 
3r     3r 

/^ »  .  • 

9r   ,  §'  _  9r 

9z'         90* 
and 

96*  _    C;'  _  36* 
r  -  -  r   -    jT 

9z  96 

With  the  above  vector  identities,  we  may  rewrite  any 

vector  relation  given  in  terms  of  x,y,z  or  r,4),z  in  terms  of 

t   «   »      t   t   I 

X  ,y  ,z   or  r  ,6  ,z  .   For  example,  it  is  easy  to  show  that 


ip'  9  ^e'  9  ^  -\  ^      ^n 

— -  r  — r  +  -r  — r  +  ez  ( — r r) 


(18) 
9r    r   96        9z 


and. 
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'"^   '^    r  9  6       90    36*  9r         9z   ae' 

+  (1,  (ilLlrl.  9t;')),-  (19) 

r     9r      96 

»   I   I 
If  we  specify  a  surface  in  r  ,0  ,z  -space  by 

«     I   t   I 

r   =  r  (0  ,z  ),  then  the  parametric  representation  of  this 

surface  is 

r((l),z)  =  r(4),z)r  +  zz 
or,  , 

»      f  t      T      »  2     '^  ' 

r(e  ,z  )  =  r^r  (0  )r   +  r^  -—  z 
%  o  o  e 

9^     .  ^^ 
Two  tangent  vectors  to  this  surface  are  — p  and  — ,  , 

9z      90 
so  that  a  normal  vector  to  the  surface  is: 

-  ^  _1 %^  ^  e_  ^  ^  j,'_  0_  9r ^  9r^  ^'  (20) 

'^        rr96     ^o9z        r98       9z 
o 

The  element  of  surface  area  is: 

I  9r     ^r   I   t   » 
dS  =  |(-^  X  ^)|d6  dz 
90     9z 

Thus,  at  the' plasma— vacuum  interface: 

dS  =  f^  d6'dz'r' \n\  [l   +   0(e^))  (21) 

It  is  easy  to  show  that  although  the  set  of  unit 

t      I      T 

vectors  x  ,y  ,z   differ  from  x,y,z,  the  corresponding  sets 
of  polar  vectors  r  ,0  ,z   and  r,(}),z  are  identical.   Thus, 
any  vector  equation  written  in  terms  of  r,$,z,  can  be  re- 
written  identically  in  terms  of  r  ,6  ,z  ,  and  for  the  re- 
mainder of  this  paper,  the  primes  will  be  omitted,  and  all 
vectors  will  be  referred  to  the  r,8,z  system  and  the  coor- 
dinates referred  to,  without  primes,  as  r,0,z.   Thus, 
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equation  (9)  can  be  rewritten  as: 

B^*^  =  YZ  +  ecr(sln  26?  +  cos  299) 

(9  ) 
-  4Aec^r^(sin  46f  +  cos  499) 

Also  (10),  (20),  and  (21)  can  be  rewritten,  dropping 

terms  of  order,  0(c)  and  0(e  ),  as: 

r  =  1  +  I  cos  26  +  (j^  -  A)c^cos  46  do') 

n  =  r  +  (c  sin  26  +  c^(|  -  4X)  sin  466)+  0(£)z     (20) 

and, 

r  , 

dS  =  j^  |n|  rd64)z  (2l) 

A  final  useful  preliminary  computation  that  follows 

from  (6),  (8),  (12),  (l4),  and  (20)  Is: 

2  2 

fi.ylO  =  ^2^2(-  c_B  (^_j,-3)  _  cB(r+r"3)  cos  29]         (22) 


where 


[B^]  =  cec6(-r^  +  r"^)  cos  26   .  (22*) 
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stability 

We  now  commence  with  the  stability  analysis  described 
in  the  Introduction,  in  which  we  study  the  second  variation 
of  the  potential  energy.   The  second  variation  of  the  poten- 
tial energy,  6W,  consists  of  three  terms,  a  volume  integral, 
6wP,  over  the  plasma  region;  a  surface  integral,  5W  ,  over 

the  plasma— vacuum  interface;  and  another  volume  integral, 

V 
6W  ,  over  the  vacuum  region.   The  system  is  considered  stable 

if  the  sum  of  the  three  contributions  is  Dositive  for  all  ad- 
missible plasma  displacements,  5(r,9,z),  and  disturbed  scalar 
potentials  of  the  vacuum  fields,  $(r,0,2),  where  $  and  C  are 

related  at  the  interface,  and  t—  =  0  on  the  bounding  outer 

on 

surface,   assumed  to  be  a  perfect  conductor. 

For  this  study,  as  we  mentioned  above,  we  are 
interested  in  determining  a  critical  beta  and  a 
corresponding,  marginally  stable  state  to  serve  as  the 
starting  point  for  bifurcation.   Thus,  we  need  only 
determine  the  lowest  eigenvalues  and  associated  eigen— 
functions  of  6W.   Therefore,  it  is  sufficient  to  determine 
£  and  $  so  as  to  minimize  the  two  volume  integrals,  holding 
n  •  C  fixed  on  the  plasma  vacuum  interface.   As  a  result,  £ 
and  $  satisfy  Euler  variational  equations,  which,  in  effect, 
reduce  the  volume  integrals  to  surface  integrals  over  the 
plasma  vacuum  interface. 

More  precisely, 

6W  =  6W^  +  6W^  +  5\i^ 
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'6WP  =  I 


dv  {  [rpCv-O^]  +  [Vx(?  X  B)f}      ,  (23) 


26W^  =  1  dS(n-5)^  n-V  ^ 
and 


=  f  dS(n-0^  n-y  if-J-      ,     (2i|) 


26W^  =  I  dv(V$)2  ,     (25) 

while  on  the  Interface: 

n-V$  =  n-Vx(?  X  B°"^)  (26) 

where  r  In  equation  (23)  Is  the  ratio  of  specific  heats. 

Provided  that  the  Euler  equations 

0    =    rpV(V-C)    -    Bx{Vx[Vx(C    X    B)]}  ,  (27) 

■AjO/a.  'x,      '\i      'x,      rx,         i\, 

0  =  V-V$  (28) 

are  satisfied,  we  may  rewrite  equations  (23)  and  (25)  as 

26WP  =   dS(n-C)[rpV-^  -  B-Vx(5  x  B)]  (23) 

26W^  =  I  dS$n'V$  (25) 

The  only  negative  term  In  6W  comes  from  equation  (24) 
[c.f.  (23)  -  (25)].   In  view  of  (22),  we  see  that  If  C  (and 

'X, 

n'C)  are  0(1)  In  e,  then  unless  all  positive  terms  In  equa— 

2 
tlons  (23)  and  (25)  are  of  order  e   or  smaller,  the  system 

Is  automatically  stable.   In  particular,  we  Infer  that 

V-C  =  0(e)  (29) 


and 


Vx(C  X  B)  =  0(e)  (30  ) 
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If  we  expand  ^  in  a  formal  series  In  e 
|(r,e,z)  =  _^Q(r,e,z)  +  e^(r,0,z)  +  £^^2^^'^'^)  +  •••   ^3l) 

then,  from  equations  (l8),  (30),  and  (32): 
9(re';(r,e,z))    3?^ 

We  may  now  rearrange  terms  among  ^   and  5,  In  equation 
(32),  leaving  C  fixed  so  that 

or  the  divergence  of  5  restricted  to  the  r,6  plane  vanishes. 
Equation  (32  )  Is  a  more  convenient  form  of  equation  (3  2). 


Euler  variational  equation,  (2  7),  as 


VJ^  =  V[pr(y-C)  -  B-Q]  =  -  Q-VB  -  B-VQ  (33) 


where.  In  view  of  equation  (3  2  ) 


Vx(^  X  B)  =  (B-V)^  -  (^•V)B  +  O(e^)      .      (3^^) 


We  set 


B  =  B   +  eB. 

A,      rv.O        'V-L 


and  find  that 

Q  =  Vx[(^^  X  B^)  +  (F   X  eB,)  +  (eC.  x  B  )  +  0{e^)-] 

» 

In  view  of  equations  (l8)  and  (32),  to  order  e, 
Q  =  ^(It  -  Ia)?.  +  e(B  -V)?^  -  (5^-V)eB^  -  Ye(y-C.)i 


20 


Denoting 
and 

then, 

-  Q-VB  -  B-VQ  =  -  e(Q^+Q, )-VB,  -  (B^+eB,  ) • V (Q^+Q,  ) 

+  0(e3)  +  o(e^)S  (35) 

Since  the  latter  expression  is  the  gradient  of  a  scalar 
function,  its  curl  must  vanish.   And,  since  the  leading 
order  term  in  this  expression  is  order  e   from  equation  (19); 
the  vanishing  of  the  r  and  6  components  of  the  curl  imply 
that  the  coefficient  of  the  z  is  of  order  e   rather  than  e  . 


Thus,  if  we  are  interested  only  in  the  leading  order  terms 

-3 

in  e,  all  terms  of  order  e   are  negligible.   Since  the  lead- 
ing order  terms  of  equation  (35)  involve  only  the  £   compo— 
nent  of  £,  we  fix  e  at  a  value,  small  enough  to  neglect 
terms  of  order  e^ ,    and  with  no  ambiguity,  refer  to  E,^   as  ? 

ryjO  % 

for  the  remainder  of  this  paper. 

We  proceed  to  expand  |  in  terms  of  the  other  small 
parameter,  c.   Since  -B-VQ  -  Q'VB  is  the  gradient  of  a 
potential  Q,,    its  curl  vanishes,  and  from  the  z-component  of 

the  curl,  we  solve  order  by  order  in  c  for  the  r,G  compo— 

r       9 
nents  C   and  5  of   E,.      These  components  are  0(1)  in  e. 
% 

Again,  from  the  vanishing  of  the  r,e  components  of 
Vx[-B'VQ  -  Q«VB],  C^  is  found  to  be  0(e)  and  hence  contrib- 


utes only  terms  of  O(e^)   to  -B'VQ 
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not  affect  this  analysis. 

More  precisely,  we  set 

5  =  £   +  cS   +  c^Cp  +  •.. 

r\j  r>^0  %1  '\,d. 

?  =  ^o  -^  ^^1  -^  ^%   ^    ■'■  (36) 

^  =  ^o  +  ^'^1  -^  ^^^2  -^  ••• 
where  the  subscripts  now  refer  to  the  order  of  c  appearing 
explicitly  in  the  expression  for  E,,    B,  and  Q. 

In  (36),  B.,  1  =  0,1,2  is  given  by  (19),  while  from  the 
definition  of  Q, 

Q.  =  I    [(B. -V)?.  .  -  (£.  . •V)B  ]  (37) 

"^   i  =  o      '^  'vJ-i     -vj-i  %  'vi 

where  B.  =  0  for  i  >  2.   From  the  Euler  equation  (33) 

{Vx  {-B-VQ-Q-VB}  }^  =  0  (38) 

Furthermore,-  since  ^  is  divergence  free,  i.e., 

'X, 

we  can  introduce  the  flux  function,  Xj  where 
^i    r  89 


w-  '       ^"^ 


2 
Xq  +  cx-L  +  c  X2  +  •  •  • 


22 


If  we  define  the  linear  operator  0^   by 

then 

O-  ^.    E  (B  .V)5 

and  the  Euler  equation  (^1)  can  be  rewritten  order  by  order 


in  c  as 


O-h'l^o  =   °  (39a) 


/T  V-Vx.  =  I    {Vx[B.-VQ.  .  +  Q.  .-VB,]} 


J  >  1     (39b) 


We  now  proceed  to  evaluate  E,.      Since  the  equilibrium, 
to  any  order  In  c.  Is  a  function  of  r  and  6  alone,  the  Euler 
equations  do  not  contain  z  explicitly,  and  we  may  thus 
assume  that  both  E,   and  i|>  have  z-dependence,  exp(-l£z). 
Unless  C   has  9  and  z  dependence,  exp[-l£(e+z) ] ,  we  see  from 
(39a)  that  ^   Is  an  harmonic  vector.   The  special  case  of  6 
and  z  dependence,  exp[-l£(9+z) ] ,  will  be  considered  separate- 
ly.  It  should  be  pointed  out,  however,  that  we  can  take 

£  7^  0  in  this  expression,  for  if  £  =  0,  regularity  of  ^^  at 

"  r  _ 

r  =  0,  and  (32)  imply  that  ^   =  0.   Thus  we  may  assume,  with 

all  generality  that  (^^,5  )  is  harmonic  in  the  plane.   Thus 

00 

E,      =      I   n&  r""^  exp[i(ne-£z)](ir  -  9) 

+      I        n5  r"^-^  exp[l(n9  -  £z)](ir  +  9)         (Ho) 

n=-l   "^'^ 

Replacing  C   in  (39b),  and  noting  that  E,      for 
n  <  0,  £  modes  are  merely  complex  conjugates  of  corresponding 
n  >  0,   £  modes,  we  have 
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/^^V-Vx,  =  l-hy^e^n{n-l){Z+n-2)6^    .r'^  ^exp  [i  f  (n-2  )  e-£z]  ]  (4l) 

Assuming, 

a    2  ^   exp[l[(n-2)e  -  iiz)  ] 

X   = J 

^  (£  +  n  -  2)"^ 

and  equating  coefficients  of  like  terms  in  (4l),  v/e  solve 

for  x-i  and  hence  5-,  .   The  result  is 

C.  =   I   n6^  „  r"^"^  exp[l((n-2)e  -  £z)][(n-2)ir  -  ne) 

^^  £  +  n  -  2 

+  (complex  conjugate  modes)  (^2) 

for  n  7^  i  1  and  £  +  n  ±  2  7^  0 .   The  arbitrary  coefficients 
6   p  will  be  determined  subsequently  so  as  to  render  the 
form  6W  stationary.   Conditions  which  render  the  Euler  equa- 
tion, (39b),  homogeneous  will  be  referred  to  as  singular  con- 
ditions, and  the  corresponding  harmonic  solutions  will  be 
called  singular  solutions.   Thus,  from  equation  (4l),  it  fol- 
lows that  n  =  i'lor£  +  ni'2  =  0  are  singular  conditions 
and  the  corresponding  solutions  are  harmonic. 

In  a  similar  way, 

/-■^V.VXo  =  I    48AY^e^n(n-l)(Jl+n-4)6^  .r"  exp  [i  ( (n-^  )  e-£z]  ] 
w  "K,  ^   '^     n=l  ' 

l-8y^e^nin-l)(n-2){i+n-h)6^    ^r'^'^exp  [i  ( (n-i|)  0-ilz)  ] 

■^=1  &  +  n  -  2 

from  which   it    follows    that    for    £+n±4/0,    £+n±2/0 

?,    =    I   Y^e^n(n-1)6       .r"^"^"  ^exp  [i  ((n-^l )  6- £z)  ]  ( (n-4  )if-ne] 

^   "^  (£+n-4)    (£+n-2) 

I   -4AY^£^n6       .r'''^^    exp  [1  ( (n-^  )  e-£z)  ]  ( (n-4  )  ir- (n+2  )  e) 


n>2 


n,£ 


£    + 


TT^ 


+  (complex  conjugate  modes)  (^3) 

Again  the  singular  solutions,  n  =  1  1,  1  2  and  £  +  n  +  4  =  0 
are  harmonic. 
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From  (19),  (^0),  (42)  and  (43),  and  in  view  of  (37), 
the  non— singular  solutions  yield 

Q^  =  I   -iYen(£  +  n)6^  ,  r"^"^  exp[i(n0  -  Jlz)](i?  -  6) 
+  (complex  conjugate  modes) 


^1  =  ° 

Q„  =   y  -iYen^6^  .r"^"^  exp[i(ne  -  ilz)](ir  -  6)  (44) 

'^^    n>2       "_12 

(£  +  n  -  2) 


+  I    iYen(n-l)6   .  r"  ^  exp[i  ((n-4) e-£z) ] ( (n-4) Gir-nS) 

n>2  ^'         

^   '^  £  +  n  -  2 

+^  -4AiYen6^  ^r"^"^^  exp[i((n-4)e  -  £z]  ]  ( (n-4  )ir-(n+2  )e) 
+  (complex  conjugate  modes) 

Next,  Q.   is  expanded  in  a  formal  series: 

2 
0      1      2 
and  computed  from  its  defining  equation,  (33).   For  the  non- 
singular  solutions,  the  result  is 

Q     =      I   iY^e^(il  +  n)^6   ,  r'^  exp[i(ne  -  £z)] 
°   n>l  "»^ 

+  (complex  conjugate) 

a.    =  l-ly^t^nia   +  n)6   „  r"^  exp[i((n-2)e  -  £z)         (45) 
^       n>l  ^'^ 

+  (complex  conjugate) 
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Q  =  I   l|XlY^E^n(ll  +  n)«„  ,  r"*2  exp[l((n-4)6  -  ilz) 

+  (complex  conjugate) 

For  the  singular  n  =  i  1  solutions 
Q   =  -  lyed   +  1)6,  „  exp[l(e  -  £z)](ir  -  9) 

^,0  1  ,  X. 

+  lyzii   -    1)6_^  ^   exp[i(-e  -  £z)(ir  +  6) 
Q,  =  -  iYeCS.  „  exp[i(-e  -  Jlz)](ir  +  6) 

+  6_^  ^  exp[i(e  -  £z)](ir  -  9)] 
Q^  =    12AiY^e^6  ,  .   r^  exp[i(3e  -  2z)](ir  -  9) 

+12AiY^e^6;L  Z     ^^    exp[i(-39  -£z)](ir  +  9) 
whereas,  for  the  singular  n  =  -  1  solutions 
Q^  =  iY^e^(£  +  1)^6-^  ^   r   exp[i(9  -  £z)] 

-  iY^e^Cii   -    D^S.i    £    I-   exp[i(-9    -    ^z)] 
Q^   =   -    2lY^e^6i    ^   r  exp[i(-e    -    S<z)] 

+    2iY^e^6_i    ^   r   exp[i(9    -    S,z)] 
^^   =    -4Xi(il-3)Y^e^<5^    _,    r^    exp[i(-39    -    Jlz)] 

-4iX(Jl+3)Y^£^5_^  _^  r^  exp[i(39  -  £z)] 
J^   =  12XiY^e^6_-^  _^  r^  expCi(9  -  Az)] 

-  12XiY^e^6;L  -1  ^^  exp[i(-e-  ^z)] 
f^^  =  i^SA^iY^e^S^  _^  r^  exp[i(9  -  ^z)] 


e+e) 


(47) 
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For  the  singular  n   =   t   2    solutions: 

Q   =  -  2±yeii    +  2)6^  „  r  exp[l(2e  -  £z)](ir  -  G) 

Q,  =  0  (^8) 

'\,1 

Q„  =  2iYe6,  .  exp[i(2e  -  «,z)](ir  -  e) 

+  2iYe62  ^   r   exp[i(-2e  -  £z)](lr  +  9)    , 

whereas,  9.    ,    Q.-,    are  non-singular  and 

Q.^   =   -ly^e^(Z   +   2)62  ^  r^  exp[i(29  -  £z)] 

-  iy^e^Cii  -  2)62  ^  r^  exp[-i(2e  +  £z)]  (^9) 

J^   =  -2iY^e^62  l^^   exp[-i£z] 

In  the  last  section  we  computed  the  equilibrium 
fields  for  both  the  plasma  and  vacuum  regions. 
In  this  section  we  have  used  these  results  together  with  the 
Euler  equations  to  evaluate  E,    Q,  and  9,.      We  now  use  these 
results  to  compute 6W^,  6W  ,  and  6W   and  finally  determine  6W 
and  its  unstable  stationary  states  to  leading  order  in  c. 

6WP: 

We  express  6wP  as  the  sum  of  three  terms 
6W  =  6WP  +   6WP  +  6W^ 
where  each  term  is  computed  as  follows: 


In  view  of  (23),  making  use  of  (^0),  (42),  (^3),    (^5), 
t         I 
(10  )  and  (20  ),  we  have  for  those  coefficients,  6  ^,    for 

which  n  7^  ±  1,  ±  2, 
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1     J     'X,   O/O     '\,l     'X,^       o      X      «; 

2  2 
=    27TY^£^H     I    (n)((£  +  n)  -  !^)  6^  ^      (50) 
n>2  ' 

while  from  (40),  (42),  and  (49),  for  6^.2  ^ 

6WP  =  2^HY^e^  J  2  ( (  £+2  )-c2)  ^6^  ^  ^^  +  2  ( (  £-2  )+c^)  26_2  ^  ^^ 


+  20^2    ^  20^2^ 


.}. 


(51) 


Finally,  from  (4o)  and  (46),  taking  into  account 
leading  order  mixed  terms  containing  the  coefficients, 

6WP  =  2TTHY^e2((Jl+l)6i^^+c6_^^^)2  +  ( (  ^"D  «_!  ^  £-^5^^  ^  J  ^ 

+  24((£-3)6_3^j^+4Ac26^^j2  +  24((5.  +  3)63^j^+4Xc26_^^j2 

(52) 


6W^ 


From  (22)  and  (H-?)^  =  {K^)'^^nl{^^  )^   =    (^l+^'\)^+0(c^) 


where  n  is  given  by  ( 20  ) ,  as  n  =  r  +  n^O  we  find  to  leading 

order  in  c     , 

2  2 

nj^o         '    n>l 

n<-l  Ji+nT'o   £,  +  n 

n>l 

3       2  1 

-J^   n^(n-2)c-^B6    5^_    ^  2c66_    6    ^       (53) 

n<-l 


The  magnetic  potential,  $,  satisfies  the  Euler  equa- 
tion (28),  and  the  boundary  condition  (26) 

at  the  plasma  vacuum  Interface  and, 

|i=  0 

dn 
at  the  outer  flux  surface. 

The  outer  boundary  of  the  vacuum  region,  i.e.,  the  con- 
ducting surface,  should  be  taken  to  be  a  flux  surface,  de- 
termined by  the  condition  i|^(r,e)  given  by  (  8  )  is  a  constant. 
However,  as  suggested  by  Weltzner  (Ref,   7 ) ,  to  simplify  the 
computations,  we  approximate  the  outer  conductor  by  the 
cylindrical  surface  r  =  p,  and  thus 

fp  (p)  =  0   .  (54) 

Since  the  boundary  condition  on  the  conducting  surface  is 
the  natural  boundary  condition  for  (28),  we  underestimate 
(overestimate)  the  stability  of  the  system  by  taking  p  so 
that  the  true  vacuum  region  lies  outside  (Inside)  the  surface 
r(0)  =  p.   In  this  manner  we  can  readily  verify  the  insen— 
sitlvlty  of  the  stability  predictions  to  the  approximate 
outer  conductor  radius,  p.   Thus,  if  we  show  that  the  system 
is  stable  for  p  less  than  some  critical  value  p  ,  the  actual 
value  of  p   can  be  taken  to  be  somewhat  larger. 
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We  proceed  with  the  computation  of  6W  .   From  (28) 


and  (5^) 


■n  ,   n  -2n. 


$  =  I    e  (r""  +  r^p"^"^)  exp[i(ne  -  £z)] 
n>o 

+  I    B  (r"  +  r-'^p^'^)  exp[i(n9  -  J^z)]  (55) 

n<o 

and, 

H-V$  =  I   n6  (-1  +  r^"p"^'^)r~""^  exp[l(ne  -  iiz)] 
n>o  '^ 

+  I  nB    a  -  r~^'^p^'')r~'^'-^   expCKnO  -  £z)]    (56) 
n<o   "^ 

From  (^0),  (42),  (43),  and  (15),  to  leading  order  In  c: 

2 
n.(Vx(?  X  B°^^))  =C£l  n6^  .r^^-^  ( (  £+n  )-  ^-]    exp[i  (ne-£z  )  ] 
'X/'V/Oi'b         n>l 

+  ((«,+l)6^  £"^^"^-1  Ji)  exp[i(e-£z)] 
+  ((£-l)6_^  ^-c&^    J  exp[l(-e  -  £z)] 
From  boundary  condition  (26),  solving  order  by  order  in  c: 

^n  =  ^%^^-  (^  '   ^))^n,£    '      ^  ^  1 

B^  =  -  ^ek-^[c6_^  j^  +  (£  +  l)<5i^j^)  (57) 

6_,  =  ?^k^(c6,^,  -  (£  -  1)6_,^J 

where , 

k   =  -g 


n    ^2n  . 

"         -       c2  -^ 
Thus,  from  (55),  (56),  and  (25),  where  n  =  (1-  77— )r  +  angular 

dependence , 

f  2 

6W^  =  27rHY^e^  |  ^  j^nL^e^  ^  J  (  £+n)-  Il|-j  2  +  l^  ( (  s,+i)  6^  ^  ^+c6_^  ^  ^]  ^ 

+  L^(-c6^^^+(£-l)6_^^j2  ^  2L2C^2^^,  2L2C^2^^^|  (58) 
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where  L^  =  k^(l  +  p"^"^). 


6W: 


24(1  +  L^)[{i  -   3)6_3^j^+  ^^c^^i^s)^ 


Combining  (50),  (51),  (52),  (53)  and  (58)  assuming  that 

2 
6  Is  of  order  c 

«W  =  ^'■HY^.^    I  {(1  .  L„)n((l  .  „)  -  ^Ul, 

n7^o,Zl,±3  •>•  ' 

+  2(1  +  L2)c^62  ^  +  2(1  +  L2)c'^6^2  ^ 

+  (l+L^)((£+l)6^^g^+c6_^^^)2  +  (l+L^)((£-l)6_^^^-c6^^^)2 

2 
^'^  ^'^   n?^o,±l   "^    "^'^ 

+  24(1  +  L2)((Jl  +  3)63^^  +  ^^c^6_-L^j^)^  I         (59) 

Note  that  for  £  +  n  =  0(1),  the  non— singular  diagonal  terms 

2 
(59)  will  dominate  unless  6   .  =  0(c  ).   Thus,  since  we  are 

interested  in  determining  the  unstable  states,  we  set 

6  .    =   0(c  )  when  i   +  n   =    0(1),  in  which  case  the  mixed 

terms  in  (53)  are  negligible  and  are  not  included  in  (59). 

It  should  also  be  noted  that  the  leading  order  compo- 
nents of  the  eigenstates  of  (59),  corresponding  to  arbitrary 
a   and  n  <  0,  can  be  identified  with  complex  conjugate  repre- 
sentations of  the  same  state  with  n  >  0.   Thus,  for 
n  ^   ±1,  ±3  we  need  only  consider  states  with  n  >  0. 

Finally,  if  the  various  components  of  the  displacement, 
C,  are  rewritten  in  the  physically  meaningful,  cylindrical 
coordinates,  the  actual  wavenumbers  i    ,    n    ,    are  given  in 
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terms  of  £  and  n  as 

I      =  —  (£  +  n) 

o 
» 
n   =  n 

If  we  relate  £*  to  the  height  of  the  plasma  column, 
i.e.,  _  2, 

X,   -  j^ 

from  (59)  each  value  of  n  corresponds  to  a  family  of  unstable 
elgenstates  under  the  conditions, 

2TTr  2Trr 

°       <  H  <  ,  n  >  2,  for  all  6 


rn   2.„  /  ng    -I        ^(-n  ^2  „  /  ng   "^ 
k   ^  ^"^(l+L^))        ^(2  ^  -V2(l+L^)-> 


2irr  2Trr  ^ 
°         <  H  <  °      ,   n=2,   6^(1+L2)c'^ 

yi+L2  VI+L2  (60) 


2^^_o ,  „  ,      ^^^o       n  =  1,   for  all  3, 

<  h.  < -       ' 


Minimizing  6W  with  respect  to  c  for  each  leading  order 
n,£  mode  of  the  above  families: 
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6W(n,£  ) 


•5  u  2  2 


[( 


Ml  +  L  ) 


1+L. 


n3£ 


)*n,. 


n  >  2,  for  all 


2TT+L^- 


53^   1 
T+lTTJ 


2,9, 


n=2,   B>(l+L2)c' 


These  instabilities  were  observed  by  Freldberg 
(Ref.   4)  in  his  study  of  helically  symmetric  systems  at 
high  betas.   They  are  referred  to  as  "kink  instabilities." 

If  one  chooses  the  column  length  to  avoid  kink  insta- 
bilities for  some  range  of  n  values  for  which  the  theory 

applies,  we  see  from  (60)  that  there  exists  a  critical  bet£ 

(L  +1)   2 
6   =  — jT c  ,  and  corresponding  to 

-(0) 


there  exists  a 


marginally  stable  state,  C 
leading  order  terms  of  X 

'-1,-1  =  • 
■^-3,-1  = 


The  coefficients  of  the 
are 


«3,-l  '  -*'=^  *l,-l 

Since  6W  is  a  homogeneous  quadratic  form,  we  can  set 
6,  _-,  =  1  in  which  case 
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X    ^°^  =  exp[i(e  +  z)](ir  -  0)  I  exp[l(-e  +  z)](i?  +  6) 

-  3Ac^r^  exp[i(-3e  +  z)](ir  +  9) 

-  SAc^r^  exp[i(3e  +  z)](lr  -  6)  (61) 

Corresponding  to  T    >  we  state  for  future  reference, 

to  leading  order  in  c, 

2 
q(0)  =  lye  c   gxp[i(e  +  z)](ir  -  9)  (62) 


and 


2 


JI^O)  =  iye  |_  r  exp[i(-9  +  z)]  (63) 

2 
■^°^  =  yek^  ^(^"^  +  ^P~^)  exp[i(e  +  z)]         (64) 
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Bifurcated  Equilibrium 

In  the  above  stability  analysis,  we  found  a  family  of 
marginally  stable  states,  E,        ,  given  by  (6l)  corresponding 
to  a  critical  beta,  3   =  (L-,+1]/^  c  .  In  this  section,  we 
will  determine  bifurcated  equilibrium  states  in  a  neighbor- 
hood of  the  original  equilibria. 

In  an  exact  bifurcation  analysis,  the  bifurcated  equi- 
libria must  satisfy  the  same  flux  constraints  as  the  original 
equilibria.   Although  we  will  not  maintain  these  constraints 
exactly,  they  will  hold  to  leading  orders  in  the  expansion 
parameters . 

To  determine  the  bifurcated  equilibria,  we  choose  one 

of  the  marginally  stable  states,  T        ,    by  arbitrarily  fixing 

the  value  of  c,  and  introducing  another  small  parameter,  d. 

We  then  expand  the  fields,  B,  plasma  surface  radius,  r,  and 

plasma  surface  normal,  n,  with  respect  to  d.   Specifically 

we  set 

B  =  b(°)  +  dB^l)  +  d^B^^)  ^  ^_^ 

'\,     O/  '\j  'V 

r  =  r(°'  -^dr<l'  +  d^r'^'  .    ...  (65) 

S=fi(°'  +dfi<^'  +d2n<2)  .  ... 

where  B    ,  r    ,  and  n    correspond  to  the  bifurcating 
state.   The  first  order  term  r    is  the  projection  of  the 
displacement,  ?    ,  onto  the  surface  normal,  n    ,  i.e., 
^(1)  ^  n^°K-^^°\      But  n^°^ 'Z^°^    is  a  complex  vector,  since 
C  has  been  expressed,  for  computational  convenience,  in 
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complex  form.   Thus,  we  may,  with  all  generality,  set 

r^l)  =  im  (n^°).l(°)) 

Thus,  from  (20)  and  (6l) 

Dos(3e  +  z) 

(66) 


r^-'-^  =  cos(e  +  z)  +  I  cos(3e  +  z) 


+  (|  -  2X)c^  cos(5e+z)  +  (^  -  A]c   cos(3e  -  z) 

Once  r    is  determined,  n    can  be  computed  from  (20) 

n^-"-^  =  sln(6  +  z)e  -  ec  sln(3e  +  z)z  (67) 

Finally,  the  first  order  perturbation  of  B,  or  first 
variation  of  B^°  ,  can  be  shown  to  be  Q^°  ,   Thus,  from  (62) 

bJJ^  =  Im  §^°'*^''  =  -  ye  ^  (sin(e  +  z)r  +  cos(e  +  z)e)(68) 

whereas,  from  (64) 

2 
g(l)out  ^    ^^  ^j  ^    ce^-  r~^[-sin(e+z)r+L,  cos(0+z)el    .(69) 

%  rv,  2  ^  1 

We  note  that  to  order  d, 

Cb2]  =  [B^°)23  +  2d[B^°)  .  B^l)] 
Since  B^  =  (L^+l)/4  c^  from  (22), 
[B^°^2]  =  i;2£2^3^V^^(_^2  ^  ^-2)  ^^^  ^9  (70) 

while  from  (9*),  (15),  (68),  and  (69): 


(o)..(l)n  -  ..2^2^3.  i  .-IrV^ 


2d[B^°^.B^-'^]  =  Y  e  c^d  <^  r-^[-i— )  cos(3e  +  z) 


(71) 


-1  ^l~l 
+  r   [-^ — )  cos(e  -  z)  +  r  cos(e  -  z)  > 
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From  (70)  and  (71),  we  verify 


2 

[B  ]  =  constant 


to  order  d  and  leading  order  In  c,  which  Is  a  necessary  con- 
dition for  the  bifurcated  fields  to  constitute  a  family  of 
equilibria. 

We  complete  the  equilibrium  analysis  by  determining  the 

2         (2)    (2)       '-  (2) 
d   terms,  B    ,  r    ,  and  n    of  the  bifurcated  expansion, 

(65).   It  Is  at  this  point  that  we  may  deviate  from  the  exact 

flux  constraints,  since  no  constraints  are  Imposed. 

(2 ) 
From  (20),  and  In  terms  of  the  unknown  r    , 

(2 ) 
n^^^  =  [-  F  +  F  cos  2(e  +  z)]?  +    [-  J     sin  2(e  +  z)  -   |f- 


36    ^ 
(72) 


In  view  of  the  above  equations  for  r,  n,  and  B,  we  can,  with 

a. 

all  generality,  set 


b[^^^   =   yecr^c,  sln(4e  +  2z) 
'v^in  1 


+  yecrcp  sin  2z  +  yecrc_  sin  26 


(73) 


while  from  (19), 

r(2)z  _ 
?ln    " 

Also,  we  can  set. 


ye  cr^c-|_  cos(ii6  +  2z)      .        (73  ) 
~2~ 


'^^^  =  a2  +  a-^  cos  2(e  +  z)  +  a^c  sin  26  sin  2(e  +  z) 
+  a|,c  cos  26  cos  2(9  +  z)  , 
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where  a,,  a^,  a^,  a^,  c,  are  determined  from  the  equilibrium 

conditions,  [b2]  =  constant  and  n  •  B  =  0. 
%  '\. 

2 
In  view  of  the  above  equations,  the  d  terms  of  n  •  B 

result  in, 

n«B  =  0  =  yecd^LC-  ^  +  a^)    sin  26 

+  (^-  +  a^  -  2a2^  +  c^)  sin  29  cos  2(6  +  z)  (7^) 
+    (2a-^^  +  2a-,  +  c-^)  cos  26  sin  2(6  +  z)] 

Next,  since  to  leading  order  in  c, 
[b2]  =  [B^°)2   ^  2d[B(°).B^l)]  +  2d2cB^°)^B^2)z   ^  ^2^^(1)2^ 

From  (70),  (71),  and  (73),  to  leading  order  in  c. 


,2,___.„„_.2  2  3^2 


(.V\.l 


[B  ]  =  constant=C  e  c  d  <  (^ — )  (i  -2a-,-  ^)    cos  26  cos  2(e  +  z) 


a- 


2  '^2      ""1   2 


L  +1      c  1 

+  (-^ — )(1+  2^)  sin  26  sin  2(6  +  z)  ^     (75) 

L  +1 
+  (-— — )(c2+c^)  sin  (e  +  z)  sin  (6-z) 

L, +1  L, -1 

+  ((-^ — )(C2-C3)  -  (-^ — ))  cos(6  +  z)cos(6-z) 

From  (68)  and  (69),  [B^^-*^]  is  found  to  be  of  smaller 
order  in  c  and,  therefore,  does  not  enter  the  present  compu- 
tation, although  it  is  stated  for  future  reference 

[B<1)2]  =  ^2^2  |1  [r-Vl  .  L?)  -  2]  (76) 

Finally,  from  (7^)  and  (75), 
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a^  =  3/4 


a2  =  a^  =  1/4 


'4  =  -  1/2 

3    2 


1   L,-l 


2 

Thus  J  we  have  succeeded  In  determining  the  d   component 

of  the  bifurcated  equilibria  which  we  summarize  as 
r^^-*  =   i  +   i  cos    2(e  +  z) 


^^   -  IT 

X-  sin  d.v    t5±n  (i^^DTz,;  -    -^ 


+  ^  sin  26  sin  2(6  +    z)    -   %   cos  26  cos  2(6  +  z) 


B^2)  =  _  2Yecr3(sin  (46  +  2z)r  +  cos  (46  +  2z)e) 


4-1.         .         ^  (77) 


+  2"  Y£cr('  -_^_y)  (sin  2zr  +  cos  2ze) 


It  should  be  noted  In  (78),  the  (sin  26r  +  cos  286) 

term  has  been  absorbed  Into  B    by  redefining  the  arbitrary 

2     ' 
constant,  c,  as  c  +  cd   =  c. 

Finally,  we  summarize  for  future  reference, 

r  =  1  +  ^  +  2.  cos  26  +  (^  -  X)c^  cos  46  +  d  cos  (8  +  z) 

+  1^  cos  (38  +  z)  +  (|  -  2X)c^d  cos  (56  +  z)       (78) 

+  (^  -  X)c^d  cos  (36  -  z) 

+  (terms  with  axial  dependence  2z) 
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From  (  78)  it  follows,  making  use  of  (20),  that 

2  2 

n  =  (1-^+1^  cos(3e  +  z)  -  1^  cos(e-z)  -  ^-   cos  2e]r 

+  ((c  -  ;^  cd^)  sin  26  +  d  sin  (6  +  z) 

^  (79) 

+  I  cd  sin  (36  +  z)  -  ^  sin  (9  -  z))9 

+  (terms  with  axial  dependence  2z) 

Finally,  making  use  of  (70),  (71),  (76),  (77)  and  (79) 
we  state  for  future  reference 

L,+3   o     L,+l 
+  [-  (^—)<i^    -    ("i— )  "]  cos  26 

3U+5 
-    ( — ^ — )  cd  cos  (6  +  z)  +  d  cos  (6  -  z)      (80) 

+  (terms  having  axial  dependence  2z)  > 
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Bifurcated  Stability 

In  order  to  study  the  stability  of  the  bifurcated  equi- 
libria, we  must  again  determine  the  sign  of  the  second  vari- 
ation of  6W  for  these  new  equilibria. 

In  the  preblfurcated  stability  analysis,  the  leading 
order  instability  was  due  to  the  singular  eigenstates, 
n  =  +1^  if,  =  -1  to  leading  order  in  c.   In  this  section,  we 
will  show  that  when  the  conducting  wall  is  sufficiently  far 
from  the  plasma  vacuum  Interface,  once  again,  the  leading 
order  instabilities  will  be  due  to  the  singular  eigenstates 

n  =  +1,  I   =  -1;  n   =    2,    I   =   -2;  and  £+n  =  ±4,  I   =    0,2  to  lead- 

2 
ing  order  in  c.   All  ^  that  are  non-singular  to  order  c   have 

been  found  to  be  stable. 

In  this  section  we  will  also  show  that  these  instabilities 
can  be  stabilized  by  placing  the  conducting  wall  within 
1.15  r   of  the  plasma  vacuum  surface. 

The  analysis  proceeds  as  before.   Stationary  values  of 
5W  are  determined  for  admissible  disturbances,  Cj  satisfying 
the  Euler— Lagrange  equations.   However,  here  we  expand  Cj  ^, 
and  Q  with  respect  to  the  parameter,  d,  in  addition  to  c. 
We  set 


Q  =  Q^°^  .   dQ^l)  .  d^Q^^)  ^ 
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%       %o       %1        0,2 
o       1        2 

O-         %0         ^1  '\,d 

for  k  =  0,1,2...  and  ^''°\    O.^^K    and  n^°''  correspond  to  the 
bifurcating  state  given  by  (6l),  (62),  and  (63),  and  as 
before  denoted  by  T        ,    Q        ,    and  ?2^°^  to  distinguish  them 
from  the  general  preblfurcated  values,  without  superscripts 

given  by  (40),  (42),  (43),  (44),  (45),  etc. 

(k) 
As   before,    Q.         Is    expressed  by    (37),    In    terms    of   the 

unknown   disturbances,    ^.       ,    and   the  bifurcated   fields    given 

'X.J 

by  (9'),  (68),  and  (77). 
If  we  set 

QUk)  ^  (  (£)    (k-£)  _r^(k-0.^.3(£) 

then  Q.    Is  expressed  as 

where, 

Bf^^=0       forl,£>2 
We  now  solve  the  Euler  equation  (38),  for  the  components 
of  £  arising  from  the  bifurcated  fields  B^  ^    and  B  K      These 
components  arise  In  two  ways:  from  the  Interaction  of  the  bi- 
furcated fields  first  with  the  bifurcating  state  X        ,    and 
then  with  the,  as  yet,  undetermined  order  d  components  of  the 
unstable  elgenstates  denoted  by  £"    .   We  proceed  with  this 

'V 

analysis . 

Due  to  the  bifurcating  state,  X        , 

Q(ii)  =  (gd)  v)r(°)  -  (I^°).y)B^i)   ,      1  =  1,2 

'\.JK       ^  ^    ^  ^  ^    ^ 
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In  view  of  (38)  and  making  use  of  (9  ),  (6l),  (62),  (68)  and 

(77),  we  find  that  due  to  the  Interaction  of  the  bifurcated 
■(o) 


U^^^'^Q^^K    Q^°).VB^^^  B^°^.Vo!,^iU  Q!,^i).VB^°)  \   =  0 
K     '^  %  %     '\.     'vijK    ,^jk    a,    J 


fields  with  ^ 

Vx 

l-x.     '^     a;     %     '\.     'vijk    ,^jk    a, 

1  =  1,2 
to  order  c   from  which  we  conclude  that  C    and  E,   ^  ^    satisfy 

the  same  equations  and  thus  have  the  same  form  up  to  order 

c  .   However,  before  we  can  conclude  that  ^^    and  ^^°''  also 

have  the  same  form,  we  must  determine  X  and  Its  Inter— 

'\, 

action  with  B    .   We  will  carry  this  out  In  the  Appendix  and 
'\. 

show  that  C^^-*  =  ^^°K 

We  determine  C"    to  leading  order  In  c  by  determining 
a. 

the  unstable  stationary  states  of  6W(C)  with  respect  to 
5  =  ^     +  dC    .   This  Is  equivalent  to  dlagonallzlng  the 


% 


% 


leading  order  terms  of  the  quadratic  form  6W(C). 

We  find  it  useful  to  decompose  6W(C)  into  the  sum  of 
four  quadratic  terms 

6w(o  =  6wcr^°\  i^°h  +  6wcr^°\  x^^h 

+  6W(C^l\  X^°h    +  6W(^^1),  C^^^) 

These  quadratics  are  further  subdivided  into  two  terms 

6w'^°'*  =  6W 
and 

6W^^^  =  6Wq^  +  5W^Q  +  5W-j_-j_ 

where  we   have    used   the   notation 

6W^.    =    6W(?^^\    ?^J^)  l,j    =   0,1  (81) 
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First  we  compute  6W    to  leading  order  in  c  and  d. 


Since 

S^ion   ._    [   r   (p^.^(o))2^.,  ip.^,^,        , 
J  Inl   %  %  '^ 

making  use  of  (80 ),  we  find  after  a  long  computation, 

6W^°)^  =  2^E^^e^l-[^^]c^-   ^  (l+L^)]  (82) 

to  leading  order  in  c  and  d. 


Noting  that 

6W^°^P  =  [  r(n.r  )r2^° 
from  (61)  and  (63) 


6W(°^P  =  2TrHc2e2^^^  3  ^4^2-,      _        ^33) 


6W(°^^: 


Again  we  consider  the  conducting  wall  to  be  a  cylinder 
of  radius  pr   in  which  case  the  vacuum  magnetic  potential 
satisfies  boundary  conditions  (26)  and  (5^)  and  is  given  to 
leading  order  by 


^^°^    =    6^(r  ^  +  rp~^)  exp[i(e  +  z)] 

+  82^^"^  +  r^p"^)  exp[iz(e  +  z)] 


(8il) 


where 


n-Q  =  n-V<I>  (26  ) 


M 


and, 

(85) 


+    Id   sln(e+z) ( 
Letting 


1   +   r^'^p    2'^)1    r""-^   expCKne    -    Jlz)] 
=   7t(o)    ,   ^(1) 


from    (  62) 


;(o)   =   , 


=    iye   —  exp[i(e    +   z)](lr   -   9)  , 


while   using    (  6l)    and    (69), 

Q^l)    =    Vx(C^°)    .    B^l)) 
'\^  a.      '^  "^ 

=  lyec   r"^  exp[12(e  +  z)](lr  +  6 ) [-^ — ) 

+    lye   |-   r    ^(L^-l)e 

Thus, 

n.V$^°)    -  2    f    .      1  3L   -1 

^  =   n-Q   =    ^ec'^    1     (-   5-  +    (— 4— )d    )    exp[l(e    +    z)] 


I   I-  2  -^   (-T 


+   ^  +   f    (1    -    2L^)    exp[21(e    +    z)] 


(86) 


From    (84),    (85),    solving    (26    )    order  by    order   In 
f  k,c2        k,c^d^  "1 

B^   =    C£  ^   4—  +  -^2 ^^2   -   l)(2k^   -   1)    j- 


62    =    Cec^d   2^   (2k^   -    1) 

where  k   and  L  were  defined  as  In  (57)  and  (58).   Since, 

6W^°^^  =  f  r|n|$^°^n.V$^°))d9dz    , 
J    %         % 

making  use  of  (84),  (85),  and  (79),  we  find  after  a  long 

computation. 
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SW'"'"    =    2i<at^z^    \    L^    =^  +   4^  [2l2(L2    -    1)    -    L^]   |>      (8?) 


6W<°' 


From    (82  ),    (83),    and    (86) 


r(o)    . 


2TTH?2e2    J    ^2^k  (_   I  +    1    [2lJ(L2    -    2)    -    L^] 


6W'^^  =  2TTH?^e^  <;  d^c   (-  ^  +  -j^  L2L^(.L2  -  2j  -  L^J  j  .-^   (88) 

The  effect  of  the  wall  on  6W   "^  can  be  determined  from 
(88)  at  once  for  the  two  limiting  cases, as  the  wall  Is 
placed  arbitrarily  far  from  the  plasma,  and  as  the  wall  is 
placed  arbitrarily  close  to  the  plasma  surface,  i.e. 
as  p  ->■  °°  and  p  ->  1. 


From  their  definitions,  lim  k  ,  lim  k^,  lim  L.^ 

p^oo         p^oo         p-voo 


and, 


from  (88  ) ,  for  p  2l  1  j 

6W^°^  =  2TrH?2^2 


(89) 


Thus  from  (88), 

lim  6W^°^  =  2^Hc^e^(-  ^  c'^d^)     .         (90) 

p^oo 

On  the  other  hand,  since 

^1 
lim  ^  =  2  ,    lim  k,  =  «> 

p-1  ^2         p-1   ^ 


L^  =  2k^  -  1    ,     L2  =  2k2  -  1     ,        (91) 


<j  (2k3  +  0(k2))c'^d2  \-  .    (92) 
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We  have  found  that  only  components  of  ^    with  coeffl— 
cients  6  ^   where  £  +  n  =  ±2,  +4,  0  and  £  =  0,  -2  result  in 
unstable  stationary  states.   The  only  6   „  with  the  above 
properties  are  6^  _^,    6^  q,  6_2  q'  *^6  -2'  "^4,0'  '^-2,-2, 
6  2i  nj  and  6    „.   We  will  show  in  Appendix  A  that  compared 
with  the  magnitude  of  the  instability  associated  with  6„  _„, 
the  other  instabilities  are  negligible.   We  will  also  show 
that  although  the  presence  of  a  conducting  wall  may  lessen 
these  instabilities,  they  cannot  be  removed.   Moreover,  we 
will  show  that  the  instability  due  to  the  leading  order  mode 
associated  with  6„   „  worsens  as  the  conducting  wall  is 
brought  closer  to  the  plasma  surface.   In  fact,  the  6^  _2 
interacts  with  ^  in  such  a  way  as  to  counteract  most  of 

the  stabilizing  effect  of  6W  ^ ^  in  the  presence  of  a  con- 
ducting wall. 

We  proceed  to  determine  6W^    due  to  the  leading  order 
eigenmode  associated  with  62  _2- 

6wii^: 

Using  the  notation  (81), 

fiW^l^P  =  6W^1)P  +  6W^1)P  +  6W^1)P 
01        10        11 


n 


where , 


01      J   ^'^  'V    ' 


6W^1^P  =  f  r(H.^^l))Ti^°)ded 


10 


a.  a- 


6W^1)P  =    r(H-5^1hf2^1^dedz 


11     J    -  ^ 

Since,  as  mentioned  above,  C^^"*,  ^^^'^  and  fi^^\  fi^°'' 
have  the  same  form,  6W   '^  is  identical  with  the  prebifurca- 
ted  value  of  5V^ .      Thus  according  to  (51), 

11  '" 

From  (63)  we  find  that 

10  ^      '" 

while  from  ( 6I)  and  (49), 

6W^^)P  =  27tHc2^2(1I  ^2^      ) 
01  "      '^'"'^ 

Thus, 

SW^l^P  =  27TH?2^2(i^^4g2^_^  +  21  ^2^^^^^^)    ^ 

(1)Z 


(93) 


6W 


From  (21)  and  (24) 


Xl)Z  _  f   r   ,-.^s2  -:„  [B^] 


6W'^'"  =    -prr  C^*^)   n-V  ^^  dSdz 
where , 

(H-O^  =  (n.(T^°^  +  ?^^h)^         .     (94) 

We  again  find  it  useful  to  decompose  6W     into 


01       10       11 
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Once  again  5W^^^^  Is  Identical  with  the  preblfurcated 
results,  and  referring  to  (53), 

6W^^^^  =  2uHc^e^( i c^sl      p)       .     (95) 

11  ^       ^'" 

From  (80)  and  (9  4),  we  find  after  a  long 

computation, 


(1)^  +  6W(1)^  =  2.H,2e2(_  li^V^  c^aS,    _,]    .    (96) 


6W 


01       10 
(1)V. 


We  decompose  6W^^   Into  three  terms: 


01        10        11 


where 


6W^1^^  =  f  r|H|$^°^n.V$^l)dedz 
01      J    '^^ 

6W^1^^  =  f  rlH|$^l^n.V$^°^dedz 
10      J    '^ 

6W^1)^  =  f  rlHI^^l^n.V^^l^dSdz 


11 

where  $    Is  given  by  (84),  while  $    is  the  scalar 
potential  that  results  from  Q    and  the  boundary  condition 
(26). 

From  (48  ), 
Q^-"-^  =  2icec^6„   „r  exp  [21(6  +  z)](ir  -  0) 

+  2iYec^62  _2i^  exp  [21(6  +  z)](ir  +  6) 
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Thus  , 


2,-2 
+   exp    [2i(-6    +    z)] 


xp    [21(6+z)] 

1 


:97) 


i 


From    (85)    and    (97),    solving    (26    )    order  by   order, 
^^^    =    Cec^62^_2   <     63(r"^   +   r^p"^)    exp    [21(9    +    z)] 

+    &^{r~^   +   r^p"^)    exp    [21(-0   +   z)] 

+    Q^ir'^   +   rp"^)    exp    [1(9    +    z)]     I 


(98) 


where. 


13   =    B^    =    cek^c    62,_2 


S(-    =    2Cec    d62   _2^^^2 


6W(1)^6W^1)^6W(1)^ 


11 


01 


10 


6W 


(1) 


Thus,    from   (8^1),    (86),    (97),    and    (98), 

r|HKn.V($^°)+$^l)))($^°U$^l))d9dz 

(99) 
2TTHc2e^(i|L262    .jc'^  +  L^  (  2L2  +  I  )c^d62   ,2) 

Adding  (93),  (95),  (96)  and  (99)  and  diagonalizlng  the 


result 


6W^^^  =  2TTH;^e2  ^  (8k2  -  k. 


\     (8k2  -  k^)c'^  [ 


2.-2 


(^(8L2-3)-H^)dn2   (^(8L2-3)-^  ^)  d^c^    "^ 
(8k2  -  k^)   J         8k2  -  k^  i 


(100) 
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As  before,  we  can  Immediately  determine  the  effect  of 
the  wall  on  the  stability  of  the  5p  _p  mode  for  the  two 
limiting  cases,  p  ^-  °°  and  p  ^  1. 

From  (100),  in  view  of  (89), 

lim  6W^^^  =  2TrH?^e^(-.246c^d^)  (101) 

On  the  other  hand,  from  (100),  in  view  of  (91),  for 
p  2:  1, 

6W^^^  =  2TTHC^e^([-  i  ^l    +    0(kj)]c^d^)  (102) 

Also,  from  ( lOQ ,  the  leading  order  part  of  the  eigen- 
vector associated  with  6„   ^  is 

K^'^''    =    26^   „r  exp  [21(6  +  z)](ir  -  9) 

where,  j- 

(^(8V3).  i^)d 
"^2,-2         8k2  -  k^ 

6W: 

As  we  previously  stated,  the  stability  of  the  bifurcated 
equilibrium  is  primarily  determined  by  the  effect  of  X  and 

the  eigenstate  associated  with  6^   _^   on  6W. 

Thus, 

where  6W*'°''  and  6W^^''  are  given  by  (88)  and  (100) 
respectively . 

It  follows  from  (91)  and  (101)  that  the  system  is  un- 
stable for  p  ^  °°j  whereas  in  view  of  (92)  and  (10  2),  the 
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system  is  stable  for  p  less  than  some  critical  value,  p^, 

where  1  <  p  <   °°  .      Table  1  shows  that  6W  _>  0  for  p  <_   I.I5, 
Thus  p   =1.15. 

p  ow  •  aw^-"-^  6w 

1.1  180  -1^2  38 

1.12  91.2  -  80.4  10.8 

1.13  71.8  -  63.6  8.2 


6W' 

:o) 

.80 

91. 

,2 

71. 

,8 

54. 

,2 

hi. 

.6 

32 

.3 

11, 

.8 

1.14  54.2         -  50.3  3.9 

1.15  41.6         -  40.8  0.8 

1.16  32.3         -  33.7  -  1.4 
1.2         11.8         -  17.2  -  5.4 

TABLE  1 
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Appendix 

We  now  verify  that,  compared  to  the  unstable  elgenmode 
associated  with  6„  _„  to  leading  order,  all  other  unstable 
states  yield  Instabilities  of  negligible  magnitudes. 
Furthermore,  In  contrast  to  the  6„  _„  mode,  the  Insta- 
bilities associated  with  these  elgenmodes  are  lessened  as 
the  conducting  wall  Is  brought  closer  to  the  plasma  surface. 

We  proceed  to  compute  6W    and  T    associated  with 
the  coefficients  6^    q,  S _^    q,  6^  q,  6_^  q,  6_^   _^,    6g  _^. 

6W^: 

Again,  since  K^-^\    5^°''  and  r/^-*,  Q^°^    have  the  same 

(1)d      '^     '^ 
form,  6W, ,  ^  Is  again  Identical  with  the  preb If urcated  value 

of  6WP  given  by  (50)  as 

6W^^^P  =  2-rTHc^e^I  inH   +   n)^6^   }  (103) 

11  n,£ 

In  Table  2,  we  list  the  leading  order  terms  of  6W_, . 
Any  term  smaller  than  leading  order  Is  considered  Insignif- 
icant for  this  analysis  and  Is  designated  by  a  cross.   We 
use  the  symbol  ■^^°''  to  refer  to  the  component  of  the  blfur- 
eating  state  with  n-th  azlmuthal  wave  number.  I.e. 
n  =  ±1,  i3. 
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COMPUTATION  OF  SW^^-'P 


-<?.-sir' 

'■<^-Li?') 

o 

2iicd 

X 

c 

4^' 

-24cd 
X 

X 

X 

«a) 

-4cd 

X 

c 

4^' 

4cd 
X 

X 
X 

«a) 

-i|cd 

X 

,(1) 

.a) 
4" 

4cd 

X 

-^2,0 

X 

X 

«a) 

(111    _    i92X)c2d 

X 

,(1) 

4" 

-48c^d 

X 

^+^,0 

4^) 

192Ac^d 

X 

"^^' 

-i|8c^d   +   c^d/2 

192Xc^d 

,(1) 

«(i) 

-^Sc^d 

X 

'^-4,0 

«") 

-192Xc^d 

X 

«<^) 

(5°5      480X)c2d 

X 

,(1) 

4" 

-2i40c^d 

X 

'>'6,-2 

4^' 

480Xc^d 

X 

"^'' 

-8c^d  -   c^d/2 

96Xc^d 

,C1) 

4^) 
4^' 

+  8c2d 

X 

^■-2,-2 

-96Xc^d 

X 

O 

0 

0 

jd) 
--2,-2 

«a) 
4^) 

0 
c^d/8    +    2c^d 

0 

X 

TABLE    2 
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Compared  to  the  leading  order  terms  In  Table  2,  we 

find  that  all  terms  from  6W   ^  are  of  lower  order.   Sum- 

10 
marizing  the  results  of  Table  2 


6W 


^)P  =  2.H,2,2  ^,  ^6^ya 


6W 


(104) 
(1)V. 


^4^ 


We  decompose  6W     as  in  (99),  and  note  that  once 

again  6W     is  identical  to  the  prebifurcated  value  given 

by  (58).   We  also  find  that  6W^^''^  +  6W^^''^  are  of  lower 

01       ^°  f  ) 
order  than  the  leading  order  terms  of  6  W   ^  .   Thus, 


6W 


6W 
(1)^. 


^^^^  =  2TTH^2e2  I  ^L^(^  +  ^)2  g2^^  |     ^     ^^^^^ 


From  (80)  and  (9^),  making  use  of  (24),  we  find  that 
6W     results  in  only  terms  of  lower  order  than  6W   ^. 


Thus, 


6W^^''^  =  0   to  order,   O(c^d^). 


6W^1) 


Adding  (10^,  (104),  and  (105),  and  diagonalizing  the 
result. 
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From  (lOQ,  the  leading  order  part  of  X  is  given  by 

^(1)  =  46^  Qr3  exp  [iiie](i?-e)  -  46_^  ^r^  exp  [-l4e](lr+e) 

-  26_2  _2^  e^P  ti(-2e  +  2z)](i?  +  6) 

+  66g  _2r5  exp  [l(6e  +  2z)](ir  -  6)  (107) 

where . 

=  _  -0568  2 
*4,0      1  +  L;^   ^  "^ 

_  _  -0039  ^2 
"^-4,0      1+Li^  ""    "^ 

*-2,-2    I+L2    ^ 

.0652  ^2, 

Also  from  (IO6)  ,  we  see  that  with  the  conducting  wall 
placed  at  a  radius  approaching  infinity,  i.e.  L^  ^  1, 

6W^^^  =  2iTHC^e2(-.3c^d2) 
Furthermore,  from  (10?)  ,  we  see  that  this  instability  is 
lessened  as  the  wall  is  brought  closer  to  the  plasma.  Refer- 
ring to  the  results  of  Table  1,  the  instability  due  to  all 
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coefficients  other  than  6p  _„  is  negligible. 

Having  computed  X        ,  we  can  now  comolete  the  computa— 

-(2) 
tion  of  5    .   As  mentioned  above,  we  must  study  the  inter— 

action  of  '^^'^^    with  b''-'-^ 

We  find  that  the  Euler  equation  yields, 
Vx  ^(B(°).V0^2))^(Q(2).,g(o))  ^  3(1).,q(1)^q(1).,3(1) 


where 


and, 


,(2)  =  (^(l).v)^(l)  _  (^(I).v)b(I) 


We  conclude  that 
C    ,  C    ,  and  ^    have  the  same  form  up  to  order  c 

%  'Xj  ^ 

-(2) 
and  therefore  the  C    component  of  the  eigenvector  does  not 

I]    2 
affect  6W  to  leading  order  0(c  d  ). 
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